Theory of nonlinear optical spectroscopy of electron spin coherence in quantum dots 
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We study in theory the generation and detection of electron spin coherence in nonlinear optical 
spectroscopy of semiconductor quantum dots doped with single electrons. In third-order differential 
transmission spectra, the inverse width of the ultra-narrow peak at degenerate pump and probe 
frequencies gives the spin relaxation time (Ti), and that of the Stoke and anti-Stoke spin resonances 
gives the effective spin dephasing time due to the inhomogeneous broadening (Ti)- The spin de- 
phasing time excluding the inhomogeneous broadening effect (T2) is measured by the inverse width 
of ultra-narrow hole-burning resonances in fifth-order differential transmission spectra. 

PACS numbers: 76.70.Hb,42.65.An, 78.67.Hc 



I. INTRODUCTION 

Electron spin coherence in semiconductor quantum 
dots (QDs) is a quantum effect to be exploited in emerg- 
ing technologies such as spin-based electronics (spintron- 
ics) and quantum computation^ The electron spin de- 
coherence is a key issue for practical application of the 
electron spin freedom and is also of fundamental inter- 
est in mesoscopic physics and in quantum physics. The 
electron spin decoherence in QDs, however, is yet poorly 
characterized. By convention, the spin decoherence is 
classified into the longitudinal and the transverse parts, 
which correspond to the spin population flip and the Zee- 
man energy fluctuation processes and are usually char- 
acterized by the relaxation time T\ and the dephasing 
time T2, respectively. Most current experiments are car- 
ried out on ensembles of spins, composed of either many 
similar QDs^^ii or many repetitions of (approximately) 
identical measurements on a single QDiLLSaLifiiiiiiSiiii 1 ^ 
The ensemble measurements are subjected to the inho- 
mogeneous broadening of the Zeeman energy which re- 
sults from the fluctuation of the QD size, shape and com- 
pound composition (and in turn the electron g-factor) 
and from the random distribution of the local Over- 
hauser field (due to the hyperfine interaction with nuclear 
spins in thermal states). The inhomogeneous broadening 
leads to an effective dephasing time T 2 * ^iLii The three 
timescales characterizing the electron spin decoherence 
can differ by orders of magnitude usually in the order 
Ti > T 2 > T 2 . For example, in a typical GaAs QD at 
a low temperature (< 4 Kelvin) and under a moderate 
to strong magnetic field (0.1 ~ 10 Tesla), the relaxation 
time T\ can be in the order of milliseconds jgjLgjjLifliil the 
dephasing time T 2 is up to several microseconds, 5,12,14 
and the effective dephasing time T 2 * can be as short as a 
few nanoseconds 

The issue is how to measure the characteristic times of 
electron spin decoherence in QDs. There have been many 
experiments both in optics 4 * 9 * 1 ^ 1 - and in transport's*^ 
which establish the spin relaxation time Ti in QDs 
of different materials. The effective dephasing time 



T 2 * is also measured for QD ensembles giving 
a lower bound of T 2 . Spin echo in microwave ESR 
experiments is a conventional approach to measuring 
the spin decoherence time T 2 excluding the inhomoge- 
neous broadening, 18,19,20 which, however, is less feasi- 
ble for III-V compound quantum dots due to the ul- 
trafast timescales in such systems (T 2 < 10 -6 sec and 
T 2 < 10~ 9 sec). Indeed, the remarkable spin echo ex- 
periments in coupled QDs done by the Marcus group are 
performed with rather long DC voltage pulses instead of 
instantaneous microwave pulses^ Alternatively, picosec- 
ond optical pulses may be used to manipulate electron 
spins via Raman processes 21 - and realize the spin echo, 
which, however, still need to overcome the difficulty of 
stabilizing and synchronizing picosecond pulses in mi- 
crosecond time-spans. A recent experiment by Greilich et 
al also shows that the inhomogeneous broadening effect 
can be filtered out from the spin coherence mode-locked 
by a periodic train of laser pulses^ 

In this paper, we will study the frequency-domain non- 
linear optical spectroscopy as another approach to mea- 
suring the electron spin decoherence times. Particularly, 
the spin dephasing rate T 2 _1 is correlated to the width of 
ultra-narrow hole-burning peaks in fifth order differential 
transmission (DT) spectra. This hole-burning measure- 
ment of the spin dephasing time is analogous to the ex- 
ploration of slow relaxation of optical coherence in atomic 
systems by the third-order hole-burning spectroscopy^ 
Here the fifth order nonlincarity is needed because the 
creation of spin coherence by Raman processes involves 
at least two orders of optical field and hole-burning two 
more. The state-of-the-art spectroscopy already has the 
ultra-high resolution (much better than MHz-resolution) 
to resolve the slow spin decoherence in microsecond or 
even millisecond timescales ^S^k 

The organization of this paper is as follow: After this 
introductory section, Sec. [H] describes the model for QD 
system and the master-equation approach to calculating 
the nonlinear optical susceptibility. Sec. IIIII presents the 
results and discussions. Sec. |W] concludes this paper. 
The solution of the master equation in frequency domain 
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is presented in Appendix 1X1 



II. MODEL AND THEORY 

The system to be studied is a semiconductor QD doped 
with a single electron. The geometry of the QD un- 
der an external magnetic field and optical excitation is 
shown in Fig. ^ ( a ) and (c). The QD is assumed of 
a shape with small thickness in the growth direction 
and relatively large radius in the lateral directions, as 
in the usual cases of fluctuation QDs and self-assembled 
QE) Sl 3 J 4 i 9 i io 4 ii i To enable the generation and manipula- 
tion of the electron spin coherence through Raman pro- 
cesses, a magnetic field is applied along a lateral direc- 
tion (x-axis). The propagation directions of the pump 
and probe laser beams are close to the growth direction 
(z-axis). The two electron spin states |±) are split by 
the magnetic field with Zeeman energy ujq. The strong 
confinement along the z-axis induces a large splitting be- 
tween the heavy hole and the light hole states, thus the 
relevant exciton states are the ground trion states |r) and 
|f) which consist of two electrons (including the doped 
one and one created by optical excitation) in the singlet 
spin state and one heavy hole in the spin state | + 3/2) 
and | — 3/2) (quantized along the z-axis with nearly zero 
Zeeman splitting), respectively. Similarly, we can also 
neglect the excitation of higher lying trions, bi-exciton 
and multi-exciton states since the energy of adding an 
exciton in each case is well separated from energy of the 
lowest trion states. The selection rules for the optical 
transitions are determined by the (approximate) conser- 
vation of the angular momentum along the growth direc- 
tion so that a circularly polarized light with polarization 
<7 + or o~ connects the two electron spin states to the 
trion state |r) or f), respectively [see Fig. Q](b)]. The 
relaxation processes in the system are parameterized by 
the exciton recombination rate Ti, the exciton dephas- 
ing rate T2, the spin relaxation rate Tj~ , and the spin 
dephasing rate . The inhomogeneous broadening 
leads to a random component e to the Zeeman splitting: 
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FIG. 1: (Color online). Schematics of (a) the quantum dot, 
(b) the selection rules for optical transitions, and (c) the op- 
tical detection geometry. 



flr(e) = e~ (eT 2 ) 2 / 2 /(\/27r/T 2 *). The hole spin relaxation is 
neglected since it is extremely slow when the hole is con- 
fined in the trion statesiii The theory presented here can 
be extended straightforwardly to include the hole spin 
relaxation, the light-hole states, the hole mixing effect 
(which leads to the imperfection in the selection rules), 
the multi-exciton states, the inhomogeneous broadening 
of the trion states, and so on, but we expect no quali- 
tative modification of the resonance features related to 
the electron spin coherence in the nonlinear optical spec- 
tra. For the interests of simplicity, we shall consider only 
cr + -polarized optical field (extension to other polariza- 
tion configurations may provide some flexibility for ex- 
periments and is trivial in the theoretical part). Thus the 
model is reduced to a A-type three-level system consisting 
of the two electron spin states |±) and the trion state |r). 
The A-type three-level model, in spite of its simplicity, 
is the basis of a wealth of physical effects including elec- 
tromagnetically induced transparency 2 ^ lasing without 
inversion^ and stimulated Raman adiabatic passage^ 
and has been successfully applied to study transient op- 
tical signals of doped quantum dots^SS 

The dynamics of the system is described in the den- 
sity matrix formalism with p a ^ being the density matrix 
elements between the states \a) and \ j3). The optical exci- 
tation and the relaxation are accounted for in the master 
equation as 



dtp T ,± = —i[£ g T ^c/2 - iT 2 )p T ,± 

-iE(t)p T , T + iE(t)p T ,± + iE(t)p±,±, (la) 
d t p T . T = -2T lPTiT + 2S[E*(t)p Tt+ +E*(t)p Tt -], (lb) 

dtP±,± = - (PtP±,± ~P±Pt-t) / T l + TlPr.r 

-2%{E*(t)p Tt± ], (lc) 
dtP±, T = ^iPt,t - i(±u c - i/T 2 )p±, T 

+iE*(t)p T ^~iE(t)p* T± , (Id) 



where £ g is the energy gap, and p± is the equilibrium 
population of the spin states in absence of the optical 
excitation, and the optical field E(t) — £\ E } e~ ia ^ 
contains different frequency components. The transition 
dipole moment is understood to be absorbed into the 
field quantities. In the rotating wave reference frame, 
the energy gap £ g is set to be zero and the optical fre- 
quencies Qj are measured from the gap. The first term in 
the righthand side of Eq. Ijldl) is the spin coherence gen- 
erated by spontaneous emission j£L2SL2L which has been 
demonstrated in time-domain experiments with signifi- 
cant effects on spin beats^ We will show that it produces 
extra resonances in fifth-order DT spectra. 

To calculate the nonlinear optical susceptibility, the 
master equation is obtained in the frequency-domain (as 
given in Appendix ^J. With the spectrum of the optical 
field given by E(Q) = 2nEj5 (Q — Qj), the density 
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derivation of the density matrix component " 
to the fifth order is lengthy but straightforward. The fina 
result is averaged with the inhomogeneous broadening 
distribution g(e). 
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III. RESULTS AND DISCUSSIONS 

The linear optical susceptibility is given by 

<?(e) 



Qj — £ g ± lu c /2 + «r 2 



de. 



(3; 



In fluctuation GaAs QDs, the exciton dephasing is mucr 
faster than the effective spin dephasing due to the inho- 
mogeneous broadening (I^ 1 < 0.1 ns < T 2 * - 10 ns)£ 
so the resonance width in linear optical spectra is usuallj 
dominated by the trion state broadening, revealing little 
information about the spin decoherence. 

In third-order optical response, the population and off- 
diagonal coherence of the electron spin are generated bj 
the Raman processes 



E 2 (2) E{ (21) . 1 • frji \-l 

P±,± * Pt'± > P±,± « (S2 2 l + 1I T \) > 



(4a) 



E 2 (2) El (21) , ;lrp\-l ( a,\ 

P±,± — > P T ,± — > P T ,± x ("21 ± + i/Ti) , (4b) 

corresponding to the illustrations in Fig. 0(a) and (b), 
respectively. Another optical field with frequency Q 1 
brings the second order spin coherence into the third or- 
der optical coherence p^±^ ■ The DT spectrum as a func- 
tion of the pump frequency f2i and the probe frequency 
f2 2 is 



(211) 
Pr,+ 



(211) 
Pt,- 



(5) 



which presents ultra-narrow resonances around f2 2 i = 
±ojo and Q 2 l — 0; with resonance width T 2 X and T-[~ , 
related to the spin population and off-diagonal coherence 
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FIG. 2: Schematics of Raman processes generating (a) the 
spin population and (b) the off-diagonal spin coherence. 



FIG. 3: (Color online). Third-order DT spectra of QDs doped 
with single electrons. The parameters are chosen such that 
the Zeeman energy u>o ~ 20 ^teV, the spin population p± = 
0.5, the pump frequency fii = £ g — 10.5 fieV, Ti — 5 peV 
(F- 1 w 0.12 ns), r 2 = 6 jueV (F^ 1 w 0.1 ns), T x = 100 ns, 
T2 = 100 ns, and T£ = 10 ns or artificially set to be 00 for the 
solid and dotted lines, respectively. The insets are enlarged 
plots showing details of the resonances. 



in Eq. (|4af> and l@bj, respectively. Such resonances are 
shown in Fig. Thus the spin dephasing time T 2 and re- 
laxation time Ti are measured. But when the inhomoge- 
neous broadening is included, since usually T2 3> T 2 * , the 
Stoke and anti-Stoke Raman resonances at 2 i — 
will be smeared to be a peak resembling the inhomoge- 
neous broadening distribution as 



( 21 ) 1 u 



-iirg (fi 2 I T w ) 



(6) 



The effect of the inhomogeneous broadening is clearly 
seen in Fig. 01 So in usual cases, the third-order DT 
spectra measure the T 2 * instead of the T 2 . The resonance 
at degenerate pump and probe frequencies (f2 2 I = 0) is 
related to the spin population and is immune to the ran- 
dom distribution of the electron Zeeman energy. So the 
spin relaxation time Ti can be deduced from the third or- 
der DT spectra, regardless of the inhomogeneous broad- 
ening. 

To measure the spin dephasing time excluding the in- 
homogeneous broadening effect, the fifth order nonlin- 
earity can be used. In the fifth order optical response, 
the spin coherence in the fourth order of optical field has 
very rich resonance structures. For instance, a double 
resonance like 



(4321 
P+,- 



(Q 3 1 -UJ c + i/T 2 ) 1 (^4321 

arises from the excitation pathway 



lo c + i/T 2 ) \(7) 



o {3) % o (3l) o (43l) 

P-,- Pt,- * P+,- Pt,- 



(4321) 
P+,- 



(8) 



as depicted in Fig. 01(a). The double resonance will mani- 
fest itself in a two-dimensional DT spectrum as an ultra- 
narrow peak at fi 3 i = ^4321 = u c with width ~ T 2 l . 
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When the inhomogeneous broadening is included, the 
ultra-narrow resonance will be smeared into a broad- 
ened peak along the direction f2 3 j = ^4321 with width 
~ I/T2 . But in the perpendicular direction (defined by 
f2 3 i = — ^432i)j the peak width remains unchanged. So 
when H4321 is fixed around u>o and f2 3 j is scanned, or vice 
versa, the DT spectrum will present a sharp peak whose 
width measures the inverse spin dephasing time T% ■ 
This peak has the character of hole-burning: The first 
frequency difference acts just as a selection of QDs with 
Zeeman energy oj c — O4321 from the inhomogeneously 
broadened ensemble. The hole-burning resonance result- 
ing from the excitation pathway in Eq. ijSJl. however, 
emerges together with the resonance associated with the 
spin population (f2 42 - + i/Ti)^ 1 as given in Eq. I|4a|) . To 
avoid the complication of mixing two types of resonance 
structures, we would rather make use another mechanism 
for spin coherence generation, namely, the spontaneous 
emission that connects the trion state to the two spin 
states through the vacuum field [related to the first term 
in the righthand side of Ea. (|TH|l ]i 3 i 29 i 30 i 31 

The generation of spin coherence in the fifth order opti- 
cal response involving the spontaneous emission can take 
a quantum pathway like 

J3) J31) J?4 (431) El (4321) r i J4321) 

P— ■ — * Ht.- "+,- Pt.- Pt.t + "—,+ ' 

(9) 

where the last step is the spontaneous emission. This 
optical process is illustrated by the Feynman diagram in 
Fig. 0] (b). The spin coherence generated by the sponta- 
neous emission and that by optical excitation can have 

opposite spin indices [p^ 1 ! — ► p^ 3 ^], which is impossi- 
ble in quantum pathways without the spontaneous emis- 
sion [as can be seen from Fig. 0] (a)]. Thus the double 
resonance becomes 

(4321) ^ IV (04321 + ^ 2r i) 

(p 3 l + lo c + iT^ 1 ) (O4321 - + iT^ 1 ) ' 

which is well separated from the spin population reso- 
nance. The spectrum is measured by fixing O4321 to be 
the hole-burning frequency uj + A with A < 1 /T 2 * and 
fine-tuning f^g to be ^4321 + S. As shown in the inset 




FIG. 4: (a) Schematics for a fourth-order optical process that 
generate spin coherence with a double resonance structure, 
(b) The Feynman diagram for the fifth order optical response 
involving the spontaneous emission, in which the optical field 
and the vacuum field are represented by the wavy arrows and 
the dotted arrow, respectively. 




A (ueV) 

FIG. 5: Contour plot of the fifth-order DT spectrum of the 
QDs as a function of A = fi 43 2i — Wo and S = Q,^ — ^4321- 
The probe frequency is fixed to be f2 4 = 22 /J,eV, the pump 
frequency Qi is fixed at 9 fieV, and the other three pump 
frequencies are scanned with f2 53 2i = (which makes = 
f2 32 = u>o + A — fi 4 j and f2 13 = Q 5 2 — ujo + A + 8, as indicated 
in the inset). The Zeeman energy u>o = 20 fieV, the spin 
population p± = 0.5, and the relaxation rates are such that 
ri = 5 ^eV (IT 1 « 0.12 ns), T 2 = 6 jiteV (rj 1 « 0.1 ns), 
Ti = 100 ns, T 2 = 100 ns, and T 2 * = 1 ns. 

of Fig. [5] the optical frequencies can be configured such 
that f^4 and Sli are fixed, fl^ are red-shifted by u>o + A + S 
from Qi, and fl§ and are red-shifted by ujq + A — f2 4 j 
from fij and O3, respectively. Thus the fifth-order op- 
tical response p^ 5 * 321 ^ oscillates at the probe frequency 
Q±, which enables the signal to be measured in the DT 
setup instead of six-wave mixing ones. We note that the 
resonance due to the spin population ~ (O5321 + i/T\) 1 
contributes only a constant background since O5321 = 
in the above frequency configuration. As shown in Fig. 
which plots the fifth order DT spectrum as a function of 
A = il 4 32i— t^o and the fine tuning 8 = ^13— ^4321, a very 
narrow hole in the spectrum as a function of f2 3 i is burnt 
around O4321, with width given by T 2 _1 . Along the di- 
rection f2 3 j = — O4321 (S = 0), the resonance is extended 
by the inhomogeneous broadening as expected. Sectioned 
plots of the DT signal with fixed A are shown in Fig. 0(a) 
for various spin dephasing time. The resonance width is 
given by the spin dephasing rate, demonstrating unam- 
biguously that the T2 is measured by the hole burning 
effect. The hole burning resonance can also be detected 
by varying the probe frequency with the pump frequen- 
cies fixed, as demonstrated in Fig. 0(b). The role of the 
spontaneous emission-generated spin coherence is veri- 
fied by the absence of the ultra-narrow resonance with 
artificial switch-off of the relevant term in Eq. (1 1 ell) . 



IV. CONCLUSIONS 

The spin coherence can be generated and detected in 
nonlinear optical spectroscopy of quantum dots doped 
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FIG. 6: (Color online) (a) The sectioned plot of Fig. with 
A = (i.e., 01321 = u) ). (b) The fifth-order DT signal 
as a function of the probe frequency with pump frequencies 
fixed to be such that fii = 9 peV, tt$ — 2 peV, and = 
Q, 5 2 — Wo = 20 ^teV. In both figures, the spin dephasing 
time T2 = 20, 50, 100 and 200 ns for the solid curves from 
top to bottom, and the dotted line is calculated with the 
spontaneously generated spin coherence artificially switched 
off (for T2 = 100 ns). The parameters are the same as in 
Fig. 



off-diagonal coherence demonstrate themselves in third- 
order differential transmission spectra as ultra-narrow 
resonances. The inhomogeneous broadening smears out 
the sharp Stoke and anti-Stoke peaks related to the off- 
diagonal spin coherence. Thus the spin relaxation time 
T\ and the effective dephasing time T 2 * are measured 
by the third-order spectra. In the fifth-order optical 
response, the generation of the spin coherence by both 
second- and fourth-order optical processes leads to dou- 
ble resonance structures in two-dimensional DT spec- 
tra, which are smeared by the inhomogeneous broadening 
along one direction in the frequency space but presents 
ultra-narrow hole-burning resonances along the perpen- 
dicular direction. So the spin dephasing time T 2 is mea- 
sured as the inverse width of the hole burning peak. 
The spontaneous emission-generated spin coherenc o 3 i 29 
is useful to produce hole burning resonances well sepa- 
rated from the spin-population resonances in the fifth- 
order spectra. The frequencies of the optical field can 
be configured properly to enable the detection of the sig- 
nal in the DT setup instead of the multi-wave mixing 
ones. In practice, the pump and probe frequencies may 
be generated from a single continuous-wave laser source 
by, e.g., acousto-optical modulation^ Since the ultra- 
narrow hole burning peaks are rather insensitive to the 
global shift of the laser frequencies and variation of the 
hole-burning frequency, non-stabilized laser sources may 
be used to resolve the slow spin decoherence^ 
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with single electrons, which is studied in this paper up 
to the fifth order nonlinearity with a A-type three-level 
model. The electron spin coherence is generated by the 
optical field through Raman processes as well as by spon- 
taneous emission of the trion. The spin population and 
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APPENDIX A: SOLUTION OF THE MASTER 
EQUATION 

The master equation in Eq. can be solved in the 
frequency domain by Fourier transformation to be 



—E(£l — ui)p± i ±(u>) — E(fl — u>)p Ti ±(u>) + E(fl — u)p TtT (oj) dto 



Q - £ g ±uj c /2 + iT 2 

+E*{n - w)pr,±(n) - E(n + wK±(n) <m 



2tT 



± J 

p±,±(u))=p±2wS(u))-(u) + iri + ip ± /T 1 ) i 

J {UJ + l/±l) (UJ + III l) Z7T 

+ (iTi - ip±/Ti) 



uj + i2Ti 2?r ' 

E*(fl - u)p T ,±(Q) - E(Q + uj)p*±(Q) dn 



UJ - UJ C + i/T 2 



(uj + i/T 1 )(uj + i2T 1 ) 
-E*(Q - w)pr,-(fl) + E(Q + oj)p*,+( n ) dtl 



uj -uj c + i/T 2 



2n 
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